NATIONAL DE 
LA PROPRIETE 
INDUSTRIE!. LE 



BREVET D' INVENTION 



CERTIFICAT D'UTILITE - CERTIFICAT D'ADDITION 



COPIE OFFICiELLE 

Le Directeur general de I'lnstitut national de la propriete 
industrielle certlfie que le document ci-annexe est la copie 
certifiee conforme d'une demande de titre de propriete 
industrielle deposee a I'lnstitut 



Fait a Paris, le . 



« 2004 



DOCUMENT DE PRIORITE 

PRESENTS OU TRANSMIS 
CONFORMEMENT A LA 
REGLE17.1.a)OUb) 



'67/141 102 



Pour le Directeur general de I'lnstitut 
national de la propriete industrielle 
Le Chef du Departement des brevets 




Martine PLANCH E 



INSTITUT 
NATIONAL DE 
LA PROPRIETE 
INDUSTRIELLE 



SIEGE 

26 Us. rue de Saint Petersburg 
75800 PARIS cedex 08 
Telephone : 33 (0)1 53 04 53 04 
Telecopie : 33 (0)1 53 04 45 23 
www.fapLfr 



BEST AVAILABLE COP\ 



ETABL1SSEMENT PUBLIC NATIONAL CR£E PAR LA LOI N« 51-444 DU 19 AVRIL 1931 



RATIO RAL DE 
LA PHOTO! ETC 

tnousTain.iv 
26 bis, rue de Saint Petersbourg 
75800 Paris Cedex 08 

Telephone : 33 (I) 53 04 53 04 Telecopie : 33 (1) 42 94 86 54 



CERTIF8CAT D'UTILITE 

Cede de la propriete inteilectuelle - Livre ^^^^ 

REQUETE EN DELI VRAWCE 
page 1/2 



N° 11354*03 



OS 540 W / 210502 



REMISE DES PIECES 

^ 20 MARS 2003 
UEU 75 !NPl PARKS 

N° O'ENREGISTREMENT Q303-4 Q 7 
NATIONAL ATTRIBut PAR UINP1 

mMttanwli 2 p MARS 2Q03 


(H NOM ET ADRESSE DU DEMANDEUR OU DU MANDATA1RE 
A QUI LA CORRESPONDANCE DOIT ETRE ADRESSEE 

Cabinet REGIMBEAU 
20, rue de Chazelles 
75847 PARES CEDEX 17 
FRANCE 

a a 


1 Vos references pour ce dossier 

\ffaadiattfj 240387 D21 012 JRC 


Confirmation d'un depot par telecopie 


□ N° attribue par IMNPI a la telecopie 


M MTURE 6E tft DEiaptff Plllfefc 




Demande de brevet * 




Demande de certificat d'utilite 


□ 


Demande divisionnaire 
1 Demande debmvet initiate 
ou demande de certijicat d'utifite' initiate 


□ 

IM° Date 1 1 1 1 1 1 1 1 1 
N° Date 1 i 1 1 1 1 1 1 1 


Transformation d'une demande de 
brevet europeen Demande de brevet initiate 


□ 

N° Date 1 I 1 I 1 I i i 1 


jgg TITRE DE L'INVENTION (200 caracteres ou espaces maximum) 



PROCEDE ET DISPOSITTF D« ATTRIBUTION DE COEFFICIENTS DE PONDERATION POUR DES CALCULS 
D'ATTITUDE 



H DECLARATION DE PRIORITE 
OU REQUETE DU BENEFICE DE 
LA DATE DE DEPOT D'UNE 
DEMANDE ANTERIEURE FRANQAISE 


Pays ou organisation 

Date| | | | | | | , | N° 
Pays ou organisation 

Date 1 i 1 i 1 t i i 1 N° 
Pays ou organisation 

Date 1 i 1 i 1 i i i 1 N° 

□ S'il y a d'autres priorites, cochez la case et utllisez Pimprime «Suite» 


gg' DEEViANDEUR (CocHez I'une des 2 cases) 


Jffl • Personne morale ' 'SM&(^S^ r}, Pe»^nne;physique^ 


Norn 

ou denomination sociale 


CENTRE NATIONAL D I ETUDES SPAHALES 


Prenoms 




Forme juridique 




N° SIREN 


1 77^S0H?l 1)1.1 


Code APE-NAF 


1 i . i 1 


Domicile 
ou 

siege 


Rue- 


2, place Maurice Quentin 

75001 PARIS / 


Code postal et ville 


1 i i i i 1 


Pays 


FRANCE 


Nationality 


Francaise 


N° de telephone (facuttati/) 


N° de telecopie (facuttati/) 


Adresse electronique (facuftatip 


□ S'il v a plus d'un demandeur, cochez la case et utilisez rimprime «Suite» 



Remplir imoerativement la 2^ jgga 




KATlOKAl 01 
LA PBOPClTtC 
IMOUSTBItlLE 



BREVET e'OSWEOT 



page 2/2 



mi 



[Reserve a HNPi > 



REMISE DES PIECES 
DATE 

UEU 20 MARS 2003 
70 W?\ PARIS 

N° O'ENREGlSTRtWEtfT 
NATIONAL ATTR1BUE PAR LIN PI 



DS 540 W/ 210502 




i N °de pouvoir permanent et/ou 
I de lien contractuel 

Rue 



Adresse 



Code postal et ville 



Pays 

N° de telephone (focuhahfj 



N° de telecopie (fiicaitaiif) 
Adresse electronique (/heuffatff) 



Les demandeurs et les inventeurs 
sont ies m§mes personnes 



20, rue de Chazelles 



7,58,47^^8 CmEKTT 



01 44-29 35 99 - 



Leiimi|^^ 



□ Oui 

® Non : Dans ce cas remplir le i ormulaire cle Designation d'inventeur(s) 



Uniau^ 



Etablissement immediat 
ou etablissement dif fere 



□ 



Paiement echelonne de la redevance 

feu denxmvemmis) 



Uniquemenl pour les personnes physiques effectual elles-snernes leur propre depot 

□ Oui 

□ Non 



H REDUCTION DU TAUX 
DES REDEVANCE8 



Uniquament pour les personnes physiques 

□ Requise pour la premiere fois pour cette invention (jowdretMaiiisdemnimpmfimO 

□ Obtenue anterieurement a ce depot pour cette invention fjoindre ime copie de /a 
decision d'admissioti a i assistance gmtuiie on indigifersa reference): AG | 1 , 



1 1 



gj SEQUENCES DE NUCLEOTIDES 
ET/OU D'ACIDES ASVltWES 



□ Cochez la case si la description contient une iiste de sequences 



Le support electronique de donnees est joint 

La declaration de conformite de la liste de 
sequences sur support papier avec le 
support electronique de donnees est jointe 



□ 
□ 



Si vous avez utilise I'imprime «Suite», 
indiquez le ttombre de pages jointes 



SIGNATURE DU DEMANDEUR 

OU DU MANDATAIRE 

(Worn et qualite du signataire) 



3 k m. 




VISA DE LA PREFECTURE 
OU DE L'INPi 



La loi ii°7S 17 du 6 janvier 1973 relative a I'informatique, aux fichiers et aux liberies s'applique aux reponses faites a ce formulaire. 
Flln jworitil nn droit d'acces el de rectification pour les donnees vous concernant aupres de I'lNPI. 



PROCEDE ET PISPO SITIF D'ATTRIBUTIQN DE COEFFICIENTS DE 
POMDERATIOM POUR PES CALCULS D'ATTITUDE. 

DOMAINE TECHNIQUE GENERAL. 

La presente invention concerne le domaine des precedes de determination 
d'orientation d'objets dans I'espace. 

Pius precisement, elle concerne un precede d'attribution de coefficients de 
masses aux mesures d'une succession d'etoiles acquises par un capteur 
stellaire afin de determiner une orientation dans I'espace. 
ETAT DE L'ART. 

De plus en plus de satellites artificiels sont en orbite autour de la Terre ou 
sont lances dans I'espace. II est bien entendu important de connaTtre la 
position dans I'espace des satellites, mais il faut connaTtre egalement leur 
attitude, c'est-a-dire leur orientation par rapport a un referentiel inertiel dans 
la voQte celeste. 

Un capteur stellaire est un equipement qui fournit sa propre mesure 
d'attitude, calculee a partir des mesures de direction d'etoiles qu'il realise, 
ou bien qui transmet directement la liste des mesures de directions- 1 
elementaires a un dispositif client, comme le logiciel central d'un satellite sur 
lequel le capteur est place par exemple. 

La figure 1 represente un satellite 1, par exemple sur une orbite 2 autour 
d'une planete, par exemple la Terre, a une altitude 13 de la surface de la 
planete. Le satellite 1 comporte sur I'une de ses parois un capteur stellaire 4 
tourne vers les etoiles de la voQte celeste, et ayant un champ de vue 
schematise par le tronc de cone 5. On a represente par le cercle 8 la partie 
de la voQte celeste observable dans le champ de vue 5. Les etoiles 10 et 11 
represented schematiquement une partie des etoiles observables par le 
capteur 4. On comprend que plusieurs, voire grand nombre d'etoiles sont 
observables en m§me temps dans le cercle 8. Les etoiles 9 et 12 
represented les etoiles de la voute celeste en dehors du champ de vue du 
capteur. 

Le satellite 1 se deplace sur I'orbrte 2 selon la fleche 6. Par consequent, on 
comprend que le cercle 8 se deplace egalement sur la voQte celeste, selon 



la fleche 7. Ainsi, Petoile 9 etait dans le cercle 8 dans le passe, et Petoile 12 
sera dans le cercle 8 dans le futur du deplacement du satellite. 
Le plus souvent, Pattitude du capteur 4 ou du satellite 1 est calculee en 
utilisant une partie seulement des etoiles presentes dans le champ de vue 8 

5 du capteur 4 stellaire. On peut par exemple choisir une dizaine d'etoiles 
pour effectuer un calcul d'attitude. Dans le calcul de Pattitude, on mesure la 
direction de Petoile 10 ou 11 par rapport a un repere 15 lie au capteur 4. 
Connaissant par ailleurs la direction de Petoile 10 ou 11 par rapport a un 
repere inertiel connu 14, on peut en deduire la direction du repere 15 - et 

10 done du satellite 1 ou du capteur 4 - par rapport au repere 14 de la voOte 
celeste. 

Le choix des etoiles a mesurer pour le calcul de la direction est determinant 
vis-a-vis de la precision de Pestimation d'attitude. 

Le choix des etoiles est important egalement quand les mesures de 
15 directions d'etoiles sont transmises directement et exploitees par le 
dispositif client du capteur stellaire. 

Dans certains cas, les etoiles ne sont pas soit rejetees soit selectionnees, 
mais une masse leur est associee, qui pondere ('importance a leur donner 
dans les traitements ulterieurs. 
20 Dans d'autres cas encore, le capteur stellaire effectue lui-meme une 
premiere selection d'etoiles, realise des mesures sur ces etoiles 
selectionnees et les transmet eventuellement accompagnees d'un poids ou 
transformees en une estimation d'attitude globale. 

Les mesures peuvent aussi avoir pour origine differents capteurs stellaires 

25 embarques a bord du meme satellite. 

Dans tous les cas, une selection est un ensemble de mesures, potentielles 
ou deja realis§es, de direction d'etoiles, mesures dont la date de validite est 
identique. Comme on Pa dit, des masses sont eventuellement associees a 
chaque mesure. Au fil du temps, un capteur stellaire 4 genere une 

30 succession de selections, tout ou partie des etoiles selectionnees faisant 
Pobjet de mesures. 

Traditionnellement, on selectionne le plus longtemps possible les etoiles 
choisies comme etant les meilleures. En effet, aux origines de la visee 
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stellaire, il etait deja difficile d'accrocher une etoile. On ne la lachait done 
plus quand on la tenait. 

Un des criteres de choix des etoiles d'une selection est notamment la 
magnitude de chaque etoile observable dans le champ de vue On choisit 
ainsi les etoiles les plus brillantes, a savoir celles qui ont la magnitude la 
plus fa.ble. Un des autres criteres de choix est la distance de I'etoile a I'axe 
optique du champ de vue. On montre sur la figure 2 plus precisement le 
cercle 8 en deplacement selon la fleche 7 sur la voute celeste. Choisir dans 
le cercle 8 des etoiles eloignees du centre optique 83 et situees dans une 
zone reference par 81 permet d'ameliorer la precision de I'estimation selon 
I'axe optique 83 du senseur. En effet, I'axe 83 est le plus sujet aux erreurs 
de mesure. Les etoiles situees dans une zone reference par 82 seront a • 
un moment ou a un autre, proches de I'axe optique 83 et provoqueront des 
erreurs de mesure. On rejette I'etoile trop proche du centre optique 83 et on 
select.onne alors une autre etoile plus eloignee du centre optique 83 
Les precedes de selection donnent actuellement une importance - 
preponderate a ces deux criteres, qui conduisent a ne pas changer 
souvent d'etoiles. 

Les procedes precedents presentent cependant des inconvenients 
En effet, certaines etoiles, dont la direction est particulierement mal estimee 
par le capteur 4, perturbent I'estimation d'attitude globale pendant toute la 
duree ou elles sent selectionnees. Cette duree peut etre assez longue 
quand I'attitude du capteur stellaire 4 evolue peu et quand son champ de 
vue 5 est plutot large. La figure 3 represente schematiquement une telle 
srtuat.on. On a trace sur la figure 3 I'erreur sur I'attitude du capteur 4 en 
fonction du temps. Les changements 30 de niveau d'erreur du graphe sont 
dus a des changements de selection d'etoiles pour la determination de 
l'att,tude. On constate que les erreurs dues au choix d'une selection sont 
relat,vement longues par rapport aux oscillations 31 sur chaque plateau de 
la courbe, qui sont quant a elles dues aux erreurs d'observation de chaque 
eto.le de la selection. On constate que les erreurs 31 s'annulent du fait de 
leur moyennage sur un temps relativement court par rapport au temps sur 
lequel est. observee chaque selection. Autrement dit, le bruit dO aux 



oscillations 31 est facilement filtrable par les techniques actuelles de 
traitement des donnees regues des capteurs stellaires par les dispositifs 
clients, car il est situe dans les hautes frequences du spectre frequentiel du 
signal acquis par le capteur. Le fait de changer peu souvent de selection 
5 provoque au contraire un bruit a faible frequence, difficilement filtrable. 

Plus generalement, les precedes d'attribution de coefficient de ponderation 
selon Tart anterieur ne permettent pas de controler la dispersion temporelle 
du bruit du a chaque selection. 

De plus, certains phenomenes - par exemple la distorsion - varient selon la 
10 position des etoiles de chaque selection dans le champ de vue. 

Plus gentralernent, les precedes d'attribution de coefficients de ponderation 
selon Tart anterieur ne permettent pas de controler ces phenomenes de 
dispersion spatio-temporelle des etoiles auxquelles sont attributes de 
masses importantes. 
1 5 PRESENTATION DE ^INVENTION. 

^invention propose de pallier ces inconvenients. 

Un des buts de I'invention est de proposer un precede d'attribution de 
coefficients de ponderation ou de masse pour un calcul d'une orientation 
dans I'espace permettant de maitriser les erreurs de mesure d'un capteur 
20 stellaire. 

Un des autres buts de I'invention est de disperser plus ou moins les etoiles 
dont les mesures de direction sont exploitees, soit par le capteur stellaire 
lui-meme, soit par le dispositif client du capteur stellaire. 
Un autre but de I'invention est de proposer un precede permettant une 

25 dispersion spatio-temporelle des selections d'etoiles. 

Enfin, un autre but de I'invention est de proposer un precede qui puisse 
prendre en consideration les selections passees afin d'etaler plus ou moins, 
sur tout le spectre, le bruit basse frequence lie aux selections d'etoiles. 
A cet effet, I'invention propose un precede d'attribution de coefficients de 

30 masses aux mesures d'une succession d'etoiles acquises par un capteur 
stellaire relie a un dispositif client afin de determiner une orientation dans 
Tespace, caracterise en ce qu'on favorise ou defavorise le renouvellement 
dans le temps, par le capteur stellaire et/ou son dispositif client, de la 
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position des mesures ayant les masses les plus importances et/ou des 
etoiles sur lesquelles portent ces mesures, afin de deplacer dans le spectre 
frequentiel une partie de !a puissance de I'erreur associee a I'ensemble des 
mesures stellaires. 

(-'invention est avantageusement completee par les caracteristiques 
su,vantes, prises seules ou en une quelconque de leur combinaison 
techniquement possible : 

- dans le calcul des masses des mesures d'une selection courante, le 
renforcement ou ('attenuation inter/ient par te biais d'un poids de distance 
associe a chaque mesure de la selection courante et caracteristique d'une 
distance moyenne entre ladite mesure d'une part et les mesures des 
selections passees et les autres mesures de la selection courante d'autre 
part ; 

- on calcule le poids de distance associe a une mesure de la selection 
courante comme une moyenne ponderee des distances respectives entre ff 
ladite mesure d'une part, et respectivement les mesures des selections': 
passees et les autres mesures de la selection courante d'autre part ; * 

- dans le coefficient de ponderation associe a la distance entre une 
premiere mesure de la selection courante et une seconde mesure d'une 
selection passee ou d'une autre mesure de la selection courante 
interviennent un coefficient de memoire associe a ladite seconde mesure, 
et/ou le poids de la seconde mesure si elle appartient a une selection 
passee ou un poids provisoire si elle appartient a la selection courante ; 

- le calcul de la distance combine la distance angulaire entre les deux 
25 mesures, et une distance identitaire fonction de la difference de nature des 

deux etoiles qui font I'objet des mesures ; 

- on definit le coefficient de memoire d'une mesure m t . a un instant t par la 
formule : 

Mem{jn t ft) = \ix ir I, - r<m '>3 , ou 
30 • T(m k ) est une date de validite d'une mesure m k . 

» u etn sont des constantes. 
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- on affecte a chaque etoile qui est I'objet d'une mesure une charge 
resumant les masses attributes aux mesures ayant porte sur ladite etoile 
par le passe, attenuees par le passage du temps ; 

- on definit la charge de Tetoile e p a un instant T par la formule : 

5 Cha(e p ,T)= j^[A{m $ )x Memoir)] 

E{m t )-e p 

0C1 Memijn^T) est le coefficient de memoire de la mesure m g a 
Tinstant T, E(mi) est 1'etoile sur laquelle portait la mesure mi et Aim^ 
est la masse de la mesure m g . 

- avant de la faire intervenir dans le calcul de la masse associee a une 
10 mesure, on actualise la charge associee a une etoile sur laquelle porte une 

mesure de la selection courante, grace a un coefficient qui depend de la 
difference A entre la date courante et la date de la derniere actualisation de 
la meme charge ; 

- que le coefficient est multiplicatif et se presente sous la forme IT* , ou n 
15 est une constante ; 

- le coefficient est additif et se presente sous la forme - P xa , ou P est une 
constante ; 

- apres avoir calculi la masse associee a une mesure de la selection 
courante, on met a jour la charge associee a Tetoile ayant fait Tobjet de 

20 cette mesure ; 

- la mise a jour est effectuee en additionnant la masse associee a la 
mesure ; 

- on fait intervenir une fonction aleatoire dans le calcul des masses ; 

- on effectue une iteration du calcul du poids de distance avec une masse 
25 provisoire pour les mesures de la selection courante, le poids de distance 

permettant de calculer une nouvelle masse qui permet elle-meme de 
calculer un nouveau poids de distance, et ainsi de suite jusqu'a 
convergence vers une masse finale ; 

- on enregistre les valeurs numeriques du procede dans des moyens 
30 memoire et de traitement du capteur et/ou du dispositif client ; 



- on augmente le renouvellement des etoiles ayant un poids important en 
elevant la frequence des mesures du capteur stellaire et/ou du dispositif 
client ; 

- on fait intervenir directement dans les masses la dispersion de I'ensemble 
5 d'une nouvelle selection grace a des moyens de traitement relies au capteur 

et/ou au dispositif client ; 

- pour faire intervenir directement la dispersion dans les masses, on utilise 
des moyens de traitement relies au capteur et/ou au dispositif client 
comportant une structure neuronale. 

1 0 L'invention concerne egalement le dispositif de mise en oeuvre du precede. 
PRESENTATION DES FIGURES 

D'autres caracteristiques, buts et avantages de l'invention ressortiront de la 
description qui suit, qui est purement illustrative et non limitative, et qui doit . 
etre lue en regard des dessins annexes sur lesquels : 
15 - la figure 1, deja commentee, represente schematiquement un satellite 
artificiel en orbite autour d'une planete et comportant un capteur stellaire ; 

- la figure 2, deja commentee, represente schematiquement un champ de' 
vue d'un capteur stellaire sur la voQte celeste ; 

- la figure 3, deja commentee, represente schematiquement un graphe de 
20 I'erreur sur I'attitude en fonction du temps dans des precedes d'attribution 

de coefficients de ponderation selon I'etat de la technique ; 

- la figure 4 represente schematiquement un graphe de I'erreur sur I'attitude 
en fonction du temps dans un precede d'attribution de coefficients de 
ponderation selon l'invention. 

25 DESCRIPTION DETAILLEE. 

Dans certaines de ses phases de calcul, un capteur 4 stellaire visible sur la 
figure 1 choisit de selectionner telle ou telle, des etoiles (10 et 11 sur la 
figure 1) qu'il detecte dans son champ de vue 8, soft pour I'utiliser lui-meme 
dans le calcul de I'attitude, soit pour transmettre directement la mesure de 

30 sa direction au dispositif client, le satellite 1 par exemple. 

La selection passe alors par ("attribution d'une masse a (Chaque etoile 10 et 
11, qui caracterise le potentiel de I'etoile vis-a-vis des besoins du dispositif 
client 1 du capteur stellaire 4. 




Selon d'autres precedes, Ies etoiles interviennent dans le calcul de I'attitude 
avec une ponderation, ou bien une masse est transmise au dispositif client 
1 avec la mesure de direction d'etoile. Toutes les etoiles sont aiors 
eventuellernent transmises, une etoile transrnise avec une masse nulle 
5 etant equivalente a une etoile non selectionnee. 

Dans chaque cas, une masse est associee, a un moment ou un autre, a 
une mesure de direction d'etoile. Cette masse caracterise le potentiel de la 
mesure vis-a-vis des besoins du dispositif client 1 du capteur stellaire 4. 

10 On va decrire un premier procede possible d'attribution de coefficients de 
masses aux mesures d'une succession d'etoiles acquises par un capteur 
stellaire afin de determiner une orientation dans Pespace. Selon ce premier 
proced§, les masses sont attributes individuellement aux mesures. 
Selon un premier procede, on controle I'importance accordee aux mesures 

15 selon qu'elles portent sur de nouvelles etoiles ou non, ce qui peut etre 
appele la dispersion temporelle des mesures d'etoiles. Elle prend en 
consideration les selections passees. On amplifie Peffet de dispersion des 
mesures. 

Selon un premier mode de mise en ceuvre possible, lors de Tobservation 
20 des etoiles, on reduit la duree consecutive pendant laquelle une meme 
etoile est selectionnee. On a represents schematiquement les effets d'une 
telle situation sur la figure 4. Ainsi, on a trace sur la figure 4 Terreur sur 
Tattitude clu capteur 4 de la figure 1 en fonction du temps. Comme sur la 
figure 3, les changements 30 de niveau d'erreur du graphe sont dus a des 
25 changements de selections d'etoiles pour la determination de I'attitude. On 
constate que les erreurs dues au choix d'une selection sont relativement 
plus courtes que ceux de la figure 3. Elles sont desormais de Tordre de 
grandeur des erreurs 31 dues aux erreurs d'observation de chaque 6toile de 
la selection. Du point de vue spectral, une frequence plus rapide de 
30 changement de selection d'etoiles revient a echanger du bruit a basse 
frequence (induit par des etoiles bruitees par les oscillations 31 et 
conservees longtemps) par du bruit a plus haute frequence (induit par un 
changement plus frequent des etoiles qui interviennent dans Testimation 



d'attitude). Plus precisement, cela revient a etaler la puissance a tres basse 
frequence dans une gamme de frequences plus large. Dans le cas oD le 
dispositif client du capteur stellaire a des capacites de mesure a haute 
frequence (ce qui est le cas quand il met en ceuvre des gyroscopes 
performants), cette operation peut etre particulierement interessante. On 
peut en effet facilement filtrer le bruit.. 

On comprend que reciproquement, si les capacites de filtrage haute 
frequence du dispositif client sont faibles, et si le bruit haute frequence est 
genant, il peut etre interessant de favoriser les basses frequences en 
conservant chaque etoile plus longtemps. On aura done un mode de mise 
en oeuvre du procede contrdlant la dispersion temporelle des mesures qui 
sera different. 

Selon un autre procede possible selon ('invention, on controle la dispersion 
spatiale des mesures d'etoiles. 

Selon un premier mode de mise en oeuvre possible, lors de ('observation 
des etoiles, on disperse les etoiles selectionnees, geometriquement dans le 
champ de vue. On contrCle ainsi la distorsion des mesures. Cette dispersion 
conceme une selection realisee a un instant donne. 

On comprend egalement que selon un autre mode de mise en oeuvre, on 
souhaite ne pas disperser spatialement les mesures. 

On combine la dispersion spatio-temporelle des mesures pour contr6ler 
I'erreur due a chaque selection dans le calcul de I'orientation. 

A cet effet, ('invention fait intervenir, dans le calcul de la masse associee a 
une mesure, la distance moyenne de la mesure correspondante avec les 
mesures des selections passees et presente. Si Ton desire que les mesures 
auxquelles a ete accordee une forte masse perdurent longtemps, on 
cherchera a minimiser cette distance moyenne. Si I'on desire au contraire 
que les mesures se renouvellent souvent on cherchera a maximiser la 
distance. 




DEFINITIONS. 

On defmit par " poids de distance " une distance moyenne d'une mesure 
avec les mesures des selections passees et presente. Le poids de distance 
interviendra ensuite dans le calcul de la masse, lequel fait intervenir 
5 d'autres considerations (par exemple la position de la mesure sur le champ 
de vue, ou la magnitude de I'etoile consideree). Une expression possible du 
poids de distance est donnee dans la suite de la presente description. On 
definit d'abord un certain nombre de ses variables possibles. 
Soft {mi,m 2 , .. m N } la suite des mesures de directions d'etoiles passees et 
10 presentes. Deux elements distincts de cette suite peuvent correspondre a la 
meme etoile physique, mesuree a deux instants differents. 
La sous-suite {mi,m 2 ,.. m P }, avec P<N, est la suite des etoiles de la 
selection presente. 

S(m i9 mj) caracterise i'ecart entre les mesures m, et rrij t et fait intervenir 
1 5 deux termes de natures tres differentes, par exemple par leur produit : 

© DA(mi,mj) - II s'agit de i'ecart angulaire dans le repere du capteur 4 
stellaire (ou tout autre distance qui s'y ramene, par exemple la distance 
geometrique sur la matrice de detection de ce meme capteur). 
Caracteriser cet ecart angulaire permet de disperser ou de regrouper, 
20 sur la matrice du capteur, les mesures auxquelles sera attribuee une 
forte masse. Cela peut permettre de s'affranchir, par exemple, d'une 
partie des erreurs de distorsion optique. 

o DE(mi,mj) - II s'agit de Tecart identitaire, c'est-a-dire de la difference de 
caracteristique des etoiles qui sont Tobjet des deux mesures. Par 

25 exemple, cet 6cart peut etre egal a 0 si les deux etoiles sont identiques, 
et a une valeur positive sinon (compte tenu du mouvement du capteur 
stellaire selon la direction 6 de la figure 1 , deux mesures correspondant 
a des etoiles identiques peuvent etre eloignees au sens de la distance 
DA). Caracteriser cette distance permet de remplacer plus ou moins 

30 souvent les etoiles qui font Tobjet des mesures. Cela permet aussi de 
favoriser des etoiles qui ont des spectres ou des magnitudes differentes 
(ou semblables) aux 6toiles precedemment selectionnees. 



Une possibility consiste a utiliser I'ecart identitaire pour caracteriser I'ecart 
entre la mesure courante et les mesures des selections passees, et I'ecart 
angulaire pour caracteriser les ecarts entre les mesures de la selection 
presente. 

5 Eventueilement, la distance peut etre une fonction de Tune ou des deux 
distances precedentes. Si la fonction est decroissante a la fois en fonction 
de DA et de DE, on calculera, au lieu d'un poids de distance, un poids de 
proximite. C'est lui qu'on pourra faire intervenir dans la formule de calcul de 
la masse au lieu ou en plus du poids de distance. 
10 Par exemple, on choisira 8{m ti m.) =qx(e D ^- m ' ) -i)x J D^(m jJ m i ) I ou q est 
une constante, DA(m,,rrtf est la distance angulaire entre les deux mesures, 
et DE(m,,mj) vaut 0 si les etoiles ayant fait I'objet des deux mesures sont 
identiques, et 1 sinon. <y(m„m y ). est done non nul seulement si les etoiles 
sont differentes. 

15 A chaque mesure m, des selections passees est associe sa masse A(m,-) 
calculee precedemment, qui fait notamment intervenir le poids de distance. ' 
Dans le calcul du poids de distance correspondant a une mesure mj de la 
selection courante, la masse des mesures passees est diminuee par le fait 
que le temps a efface la memoire que Ton garde de cette mesure. Le 

20 resultat de cette diminution est la masse remanente. 

La masse remanente peut etre calculee comme le produit de la masse et 
d'un coefficient de memoire Mem(p h /t) qui caracterise la memoire que Ton 
garde, de la mesure m, a la date t. A la date T( mj ) de la mesure mj , le 
coefficient de memoire de la mesure m ( s'ecrit Mem{pi t IT{tn^ . 

25 Le coefficient de memoire Mem(m. / Odecroit quand la mesure m, s'eloigne 

dans le passe. Cela traduit I'effacement progressif de la memoire des 
mesures stellaires. 

On retiendra par exemple le produit suivant : 
Mem(m i lt) = ]xx n~ [t ~ r( ' n ' )J ; oil 
30 • T(m k ) est la date de validite de la mesure m k . ' 
® u et n sont des constantes. 
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Les valeurs des constantes peuvent dependre de I'ecart S(m i9 mj) . Ainsi, en 
faisant varier p, le coefficient de nnemoire peut passer a la valeur 0 quand 
[t-Tinii)] depasse une certaine valeur. Cela efface la memoire de la 
mesure. Le coefficient peut aussi passer a la valeur 0 quand [t - Hm,)] est 

5 inferieur a une certaine valeur. Cela permet de forcer la selection d'etoiles 
presentes dans les selections tres recentes. Cela permet aussi d'exclure de 
la formule du poids de distance les mesures de la selection courante. 

CALCUL DU POIDS DE DISTANCE. 

Le poids de distance associe a la mesure mj est une moyenne, sur 
10 Tensemble des mesures des selections passees et des mesures de la 
selection courante (sauf mj). 

La moyenne peut §tre ponderee par les masses remanentes des mesures 
qui interviennent dans la moyenne. 

La moyenne sera calculee par exemple sous la forme : 

K 



15 P{mj) = 



20 



M 



; ou 



€1 et q> sont des constantes. Si Q et <p valent 1 il s'agit d'une moyenne 
lineaire des distances. Si 9 vaut 2 il s'agit d'une moyenne quadratique de 
ces memes distances (Q vaudra en principe 1). Si <p<l, on accentue 
Teloignement qui concerne les paires de mesures tres proches. 



Si la mesure mj fait partie de la selection presente, la masse A(mj) ne fait 
pas intervenir le poids de distance. 

Dans la pratique, la masse d'une mesure de la selection courante peut etre 
calcutee & partir d'un poids de distance arbitraire. Suite a ce premier calcul, 
25 une methode iterative peut etre mise en ceuvre : a chaque fois intervient, 
dans le calcul des poids de distance des mesures de la selection courante, 
les masses respectivement associees a ces mesures, et qui prennent en 
compte le poids de distance calcute a Titeration precedente. Tous les poids 
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de distance des mesures de la selection courante sont alnsi recafcules 
jusqu'a leur convergence. 

Une autre possibilite consiste, dans le calcul des poids de distance de la 
selection courante. a ne pas faire intervenir les autres mesures de la 
5 selection courante. Cela revient a choisir Mem( mi /t)=0 pour toutes les 
valeurs de i telles que mj soit une mesure de la selection courante. 

Le calcul du poids de distance permet d'obtenir une masse pour chaque 
mesure d'etoile. 

1 0 Par exemple, on calculera la masse comme A{ mj ) = e M x ^(« y )x?( m .) , ou 
M est la magnitude .de I'etoile qui fait I'objet de la mesure mj, v ( mj ) la 
distance angulaire entre la mesure mj et le J'axe optique du capteur stellaire, 
et P(mj) le poids de distance associe a la mesure m L 

On favorise ainsi dans une selection les etoiles ayant les masses les plus ' 
15 importantes. 

On comprend alors que I'on controle la dispersion spatio-temporelle des 
selections en jouant sur les differentes variables du poids de distances, qui 
donne une masse finale differente aux mesures. On change plus ou m'oins 
vite les selections, et les selections sont plus ou moins dispersees entre 
20 elles sur la matrice du capteur, en fonction des parametres des variables 
que I'on choisit. 

On va maintenant decrire une variante du procede selon lequel, pour limiter 
('information a memoriser et pour simplifier les calculs, le poids de distance 
25 des mesures des selections passees est calcule recursivement. 

Pour cette variante du procede, on utilise la distance identitaire pour definir 
la proximite entre les mesures. 



DEFINITIONS. 

30 Selon la variante du procede, I'information qui resume le passe caracterise 
non pas directement les mesures elles-memes, mais les etoiles sur 
lesquelles ont porte les mesures des selections passees. Cette information 



se resume a un couple C p assode a chaque etoile e p dont la mesure a fait 
partie d'une selection passee d6fini par : 
C p ={Cha,date}, ou 

date est la date de validite de la charge Cha, laquelle evoiue dans le temps. 
5 La charge Cha est le cumul, a I'instant considere, des masses remanentes 
associes a ('ensemble des mesures qui, dans les selections passees, ont 
porte sur I'etoile e p . La charge associee a I'etoile e p a un instant T se definit 
done theoriquement, selon les notations definies plus haut et si E(nrii) est 

I'etoile sur laquelle portait la mesure m i( par : 

// 

10 Cha{e p ,T) = ]T [Aim,) x Mem^ I T)\ 

Pour prendre en consideration la dispersion spatiale entre les selections 
successives (la zone visee dans le del pouvant evoluer sous Teffet du 
mouvement du satellite), on pent multiplier le terme a Tinterieur de la 
somme par une fonction de la distance angulaire entre la mesure courante 
1 5 portant sur I'etoile e p , et la mesure m { . Cette fonction prendra par exemple la 
forme : 

-i/xdMep),™,) 

ou m(e p ) est la mesure de la selection courante portant sur I'etoile p, d la 
distance entre les deux mesures en argument, et tj un coefficient. 
20 A la place d'une charge, on peut calculer un souvenir associe a Tetoile, 
moyenne des coefficients de memoire associes aux mesures ayant porte 
sur Tetoile consideree, ponderes par leurs masses, par exemple : 

^[AimJxMemim, IT)} 
SOU{e p ,T) = ^ B 

25 

Ce souvenir est une « charge normalisee » qu'on fera intervenir dans le 
calcul de la masse. 



Propagation des charges. 

Au moment du traitement des mesures de la selection courante, on propage 
les charges de la facon suivante. 

On effectue un rafraTchissement des charges et une mise a jour des dates 
5 associees. Les charges associees aux etoiles qui font partie des selections 
passees sont diminuees pour traduire I'effacement de la memoire. 
Eventuellement, cette mise a jour peut concemer seulement les etoiles de 
la selection courante, mais les informations concernant les autres etoiles 
doivent toutes etre conservees. Les charges associees aux etoiles ne 
1 0 faisant partie d'aucune selection passee restent nulles. Les dates associees 
aux charges des etoiles des selections passees sont passees a la date de 
la selection courante. 

Cette diminution de la charge peut s'effectuer, si A est I'ecart temporel defini 
comme suit : 

15 A= date de la selection courante- date precedemment associee a la 7 

charge, 

par Tun des deux procedes suivants : 
1. multiplication de la charge precedente par un facteur : 

ET A , ou nest une constante 
20 2. addition a la charge precedente d'un facteur : 

- P xa , ou P est une constante. 
La charge ainsi obtenue est forcee a 0 si le resultat de I'operation 
precedente est negatif. 

La charge est calculee sans approximation, quand a la fois on a choisi le 
25 premier precede et que le coefficient de memoire se presente sous forme 
exponentielle. Dans les autres cas, le calcul est approche. 
Ensuite, on attribue des poids de distance. On attribue comme poids de 
distance a chaque mesure de la selection courante, la charge rafraTchie 
associee a I'etoile qui fait I'objet de la mesure consideree. Si I'etoile n'avait 
30 fait I'objet d'aucune mesure precedente, sa charge de' depart est nulle et 
son poids de distance le devient done aussi. 
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On calcule ensuite les masses. Les poids de distance calcules lors de 
I'etape precedente permettent le calcul des masses associees aux 
mesures de la selection courante, par exemple grace a la formule 
A(nij) = e~ M x y/{mj)xP(mj) . 

5 On actualise ensuite les charges. Les charges des etoiles sur lesquelles 
portent les mesures de la selection courante sont augmentees de la valeur 
des masses respectives des mesures. Pour chaque etoile e h qui n'etait pas 
precedemment caracterisee par un couple C h ={Cha,date}, la date est fixee 
a la date de la selection courante. Pour toutes les etoiles concemees, 
10 I'affectation de la date peut etre realisee a cette etape plutot qu'a I'etape 
rafraTchissement des charges. 

Enfin, on annule les charges trap faibles. Cette etape est optionnelle. Les 
charges qui sont inferieures a une valeur plancher peuvent etre annulee, ce 
qui revient a oublier que I'etoile a fait I'objet d'au moins une selection par le 
15 passe. 

Dans le calcul des masses, on peut egalement faire intervenir une fonction 
aleatoire, par exemple une variable aleatoire gaussienne ou uniforme. Cela 
contribue a disperser spatio-temporellement les etoiles, sous I'effet du 
hasard cette fois. 

20 MlSE EN CEUVRE DU PROCEDE. 

Le calcul du poids de distance ou du couple de charge suppose de 
memoriser des informations relatives aux mesures des selections 
anterieures, choisies par exemple, pour chaque mesure, parmi : 



© La date de la mesure, 



25 



© La masse finalement attribute a la mesure, 



La direction de la mesure, reperee angulairement dans un repere 
15 lie au capteur 4 ou vectoriellement dans un repere lie a sa 
matrice, 



Certaines caracteristiques de I'etoile qui a donne lieu a la mesure. 




Dans le cas de I'utilisation du procede recursif, on garde le couple 
C p ={Cha,date} pour chaque etoile ayant fait I'objet d'en ensemble de 
mesures dont la memoire n'est pas effac^e. 

Des moyens memoire et de traitement sont done associes au capteur et/ou 
5 a son dispositif client, a savoir le satellite. 

Dans le cas oD la limite de la memoire disponible est atteinte, on peut 
supprimer les informations concernantles mesures dont la memoire est la 
plus faible. Le critere pourra porter sur : 

* le produit de la m§moire associee a la mesure par la masse de 
10 cette meme mesure. 

© la charge, en particulier en cas d'utilisation du procede recursif. 

On va decrire un deuxieme procede possible d'attribution de coefficients de ; 
masses aux mesures d'une succession d'etoiles acquises par un capteur^ 

15 stellaire afin de determiner une orientation dans I'espace. Selon ce£ 
deuxieme procede, on choisit directement la meilleure selection d'6toiles 
pour la determination de Torientation dans I'espace. On ne considere plus 
les etoiles individuellement pour determiner celles qui ont les plus gros 
coefficients, avant des les regrouper dans des selections. On fait intervenir 

20 directement dans les masses la dispersion de Tensemble d'une nouvelle 
selection grace a des moyens de traitement relies au capteur. Selon ce 
procede, les moyens de traitement relies au capteur peuvent comporter une 
structure neuronale. La structure neuronale permet d'attribuer des masses 
aux etoiles pour la determination de I'orientation du capteur ou du satellite, 

25 apres un processus d'apprentissage par exemple. Le critere 
d'apprentissage du r6seau neuronal peut faire intervenir les coefficients de 
distances, distances moyennes, charges et/ou souvenir d§finis 
precedemment. II peut aussi faire intervenir la moyenne, sur les etoiles 
d'une selection, des poids de distance eventuellement ponderes par les 

30 masses des etoiles (ou par les masses provisoires qui ne tiennent pas 
compte des poids de distance). 




Dans les deux precedes, le phenomene de dispersion des etoiies 
selectionnees peut etre amplffie en forgant ie capteur steliaire a travailler a 
une frequence plus elevee que la frequence strictement utile. En effet, 
moyennant la mise en osuvre des proGedes precedents le nombre de 
5 selections utilisees sera augments, d'ou un etalement supplemental du 
bruit sur le spectre. Les avantages ainsi obtenus sont de nature a 
compenser, du point de vue du dispositif client du capteur steliaire, 
I'augmentation de bruit induit par le raccourcissement de Integration 
d'information lumineuse par le capteur steliaire. 
10 Avant utilisation, les mesures stellaires peuvent alors etre filtrees puis sous- 
echantillonnees, pour revenir a la frequence souhaitee. Le filtrage passe en 
principe par le calcul du quaternion d'attitude. La methods est efficace 
meme si Ton n'utilise pas, dans le critere, les poids du passe et de 
proximite. 

15 

Les developpements qui precedent s'appliquent aussi bien d'une part a un 
ensemble d'etoiles detecte et selectionne' regulierement (par exemple une 
fois par seconde), au meme moment - comme e'est le cas pour un 
fonctionnement classique du capteur steliaire a detecteur CCD -, ou d'autre 
20 part a ensemble d'etoiles ou les etoiies sont detectees irregulierement 
(utilisation de nouveaux types de detecteurs) et oD les selections sont 
produites irregulierement, comme e'est le cas des capteurs stellaires a APS 
(Active Pixel Sensor). 

Les precedes peuvent etre mis en oeuvre a plusieurs niveaux, 
25 eventuellement simulianement, par exemple : le niveau de la selection des 
etoiies sur lesquelles le ou les Gapteurs stellaires vont faire porter leurs 
mesures ; la ponderation des mesures effectuees au niveau du ou des 
capteurs stellaires ; la ponderation des mesures effectuees au niveau du 
dispositif client des mesures. 

30 



REVINDICATIONS. 

1. Precede d'attribution de coefficients de masses aux mesures d'une 
succession d'etoiles acquises par un capteur stellaire (4) relie a un 
dispositif (1) client afin de determiner une orientation dans I'espace, 
caracterise en ce qu'on fayorise ou defavorise le renouvellement 
dans le temps, par le capteur (4) stellaire et/ou son dispositif (1) 
client, de la position des mesures ayant les masses les plus 
importantes et/ou des etoiles sur lesquelles portent ces mesures, afin 
de deplacer dans le spectre frequentiel une partie de la puissance de 
I'erreur associee a I'ensemble des mesures stellaires. 

2. Precede selon la revendication 1, caracterise en ce que, dans le 
calcul des masses des mesures d'une selection courante, le 
renforcement ou I'attenuation intervient par le biais d'un poids de 
distance associe a chaque mesure de la selection courante et 
caracteristique d'une distance moyenne entre ladite mesure d'une 
part et les mesures des selections passees et les autres mesures de 
la selection courante d'autre part. 

3. Precede selon la revendication 2, caracterise en ce qu'on calcule le 
poids de distance associe a une mesure de la selection courante 
comme une moyenne ponderee des distances respectives entre 
ladite mesure d'une part, et respectivement les mesures des 
selections passees et les autres mesures de la selection courante 
d'autre part. 

4. Precede selon la revendication 3, caracterise en ce que, dans le 
coefficient de ponderation associe a la distance entre une premiere 
mesure de la selection courante et une seconde mesure d'une 
selection passee ou d'une autre mesure de la, selection courante, 
interviennent un coefficient de memoire associe a ladite seconde 
mesure, et/ou le poids de la seconde mesure si elle appartient a une 




selection passee ou un poids provisoire si elle appartient a la 
selection courante. 

5. Procede selon Tune des revendications 3 ou 4, caracterise en ce que 
le calcul de la distance combine la distance angulaire entre les deux 
mesures, et une distance identitaire fonction de la difference de 
nature des deux etoiles qui font I'objet des mesures. 

6. Procede selon I'une des revendications 4 ou 5, caract6rise en ce 
qu'on d6finit le coefficient de memoire d'une mesure m t a un instant 
t par la formule : 

MemOti; /t)=nx n~ [,-r(mi) l , ou 
e T(m k ) est une date de validite d'une mesure m k . 
o p et n sont des constantes. 

7. Procede selon I'une des revendications 2 a 6, caracterise en ce 
qu'on affede a chaque etoile qui est I'objet d'une mesure une charge 
resumant les masses attributes aux mesures ayant porte sur ladite 
ttoile par le pass<§, attenutes par le passage du temps. 

8. Selon la revendication 7, caracteris6 en ce qu'on definit la charge de 
I'etoile e p a un instant T par la formule : 

Cha{e p ,T)= Y\A{m,)x Memim^T)] 

oil Memitn.JT) est le coefficient de memoire de la mesure m, a 
I'instant T, E(mi) est l'6toile sur laquelle portait la mesure mi et A(m,) 
est la masse de la mesure . 

9. Procede selon I'une des revendications 7 ou 8, caracterise en ce 
que, avant de la faire intervenir dans le calcul de la masse associee 
a une mesure, on actualise la charge associee a une etoile sur 
laquelle porte une mesure de la selection courante, grace a un 



coefficient qui depend de la difference A entre la date courante et la 
date de la derniere actualisation de la meme charge. 

10. Precede selon (a revendication 9, caracterise en ce que le coefficient 
est multiplicatif et se presente sous la forme n~ A , ou IT est une 
constante. 

1 1. Precede selon la revendication 9, caracterise en ce que le coefficient 
est additif et se presente sous la forme - px &, ou P est une 
constante. 

12. Precede selon I'une des revendications 7 a 11, caracterise en ce 
que, apres avoir calcule la masse associee a une mesure de la 
selection courante, on met a jour la charge associee a I'etoile ayant 
fait I'objet de cette mesure. 

13. Precede selon la revendication 12, caracterise en ce que la mise a 
jour est effectuee en additionnant la masse associee a la mesure. 

14. Precede selon I'une des revendications 1 a 13, caracterise en ce 
qu'on fait intervenir une fonction aleatoire dans le calcul des masses. 

15. Precede selon I'une des revendications 1 a 14, caracterise en ce 
qu'on effectue une iteration du calcul du poids de distance avec une 
masse provisoire pour les mesures de la selection courante, le poids 
de distance permettant de calculer une nouvelle masse qui permet 
elle-meme de calculer un nouveau poids de distance, et ainsi de 
suite jusqu'a convergence vers une masse finale. 

16. Precede selon I'une des revendications 1 a 15, caracterise en ce 
qu'on enregistre les valeurs numeriques du precede dans des 




rnoyens memoire et de traitement du capteur (4) et/ou du dispositif 
(1) client. 

17. Precede selon Tune des revendications 1 a 16, caracterise en ce 
5 qu'on augmente le renouveilement des etoiles ayant un poids 

important en elevant la frequence des mesures du capteur stellaire 
(4) et/ou du dispositif (1) client. 

18. Procede selon la revendication 1, caracterise en ce qu'on fait 
10 intervenir directement dans les masses la dispersion de Tensemble 

d'une nouvelle selection grace a des rnoyens de traitement relies au 
capteur (4) ettou au dispositif (1 ) client. 

19. Procede selon la revendication 18, caracterise en ce que pour faire 
15 intervenir directement la dispersion dans les masses, on utilise des 

rnoyens de traitement relies au capteur (4) et/ou au dispositif (1) 
client comportant une structure neuronale. 

20.Systeme de poursuite ou d'acquisition d'etoiles comportant un 
20 capteur stellaire (4) relie a un dispositif (1) client comportant des 

rnoyens d'attribution de coefficients de masses aux mesures d'une 
succession d'etoiies acquises par le capteur ou son dispositif client 
afin de determiner une orientation dans Pespace, caracterise en ce 
qu'il comporte des rnoyens aptes a favoriser ou defavoriser le 
25 renouveilement dans le temps, par le capteur stellaire (4), des 

mesures ayant les masses les plus importantes, et des rnoyens 
aptes a deplacer dans un spectre frequentiel une partie de la 
puissance de Terreur associee a I'ensemble des mesures stellaires. 

30 21.Systeme selon la revendication 21, caracterise en ce qu'il comporte 
des rnoyens relies au capteur (4) et/ou au dispositif (1) client 
comportant une structure neuronale. 
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